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Abstract 

The asymmetric simple exclusion process (ASEP) on a one-dimensional lattice is 
a system of particles which jump at rates p and 1 — p (here p > 1/2) to adjacent 
empty sites on their right and left respectively. The system is described on suitable 
macroscopic spatial and temporal scales by the inviscid Burgers' equation; the latter 
has shock solutions with a discontinuous jump from left density p_ to right density 
p + , p_ < p + , which travel with velocity (2p — 1)(1 — p + — />_). In the microscopic 
system we may track the shock position by introducing a second class particle, which 
is attracted to and travels with the shock. In this paper we obtain the time invariant 
measure for this shock solution in the ASEP, as seen from such a particle. The mean 
density at lattice site n, measured from this particle, approaches p± at an exponential 
rate as n — > ±00, with a characteristic length which becomes independent of p when 
p/(l—p) > yj p + {l — /9_)//9_(l — p + ). For a special value of the asymmetry, given by 
p/(l — p) = — ~~ P+)i the measure is Bernoulli, with density p__ on the 

left and p+ on the right. In the weakly asymmetric limit, 2p — 1 — ► 0, the microscopic 
width of the shock diverges as (2p — 1) . The stationary measure is then essentially 
a superposition of Bernoulli measures, corresponding to a convolution of a density 
profile described by the viscous Burgers equation with a well-defined distribution for 
the location of the second class particle. 
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1. Introduction 



The aim of this and our previous work [[T], |2], with S. Janowsky, is the determination 
of the underlying microscopic structure of a fluid in regions in which it has shocks on the 
macroscopic scale. We consider situations in which the system evolves macroscopically 
according to some deterministic autonomous equations, such as the Euler, Navier-Stokes, 
or Burgers' equations ||, |J |5[]. When the solutions of these equations are smooth, we 
can assume that on the microscopic level the system is (essentially) in a local equilibrium 
state, determined by the local macroscopic parameters obtained from the solutions. What 
is less clear, however, and is of particular interest, is what happens when the macroscopic 
evolution is not smooth, as in the occurrence of shocks; these are described by regions 
of very large gradient in solutions of the Navier-Stokes or viscous Burgers' equations and 
by discontinuities in solutions of their zero viscosity limits, the Euler or inviscid Burgers' 
equations. It is first of all not clear in what sense the macroscopic equations are to be 
interpreted in such regions, since their derivation (heuristic or rigorous) is based on the 
assumption of slow variation in the system's properties on the microscopic scale. Beyond 
that, these equations do not describe the structure of the shocks on the microscopic scale. 
This is the problem we wish to address here. In particular, do the statistical properties 
of the atoms or molecules change abruptly on the interparticle distance scale? Or do 
these properties change significantly only over much larger distances? (For a discussion of 
shock structure in gases based on the "mesoscopic" description provided by the Boltzmann 
equations, see 0.) 

One model for which these questions have been answered, at least partially, is the 
asymmetric simple exclusion process (ASEP) || |J on the one- dimensional lattice In 
this model particles attempt at random times (distributed as independent Poisson processes 
of unit density at each site) to jump to an adjacent site, choosing the site on their right 
with some fixed probability p and that on their left with probability q = 1 — p; the attempt 
succeeds if the target site is not already occupied. For any value of p the set of extremal 
translation invariant stationary states is the set of Bernoulli measures v p (0 < p < 1) [H]; 
in the state z/ p , each lattice site is occupied independently with probability p and there is 
therefore a current {2p — l)p(l — p). The dynamics satisfies detailed balance with respect 
to z/p if and only if the transitions are symmetric (p = 1/2). 

The macroscopic mass density u(y,t), which is an appropriately scaled continuum 
limit of the particle density in the microscopic model, is described PH-fll^l by the inviscid 
Burgers' equation 





1 



It is well known || that solutions of ( |1 . 1| ) can exhibit shocks. In the simplest example, 
u(y,t) = p + for y > y^it) and u(y,t) = p_ for y < yo(t), where < p_ < p+ < 1; the 
shock position yo(t) moves with the constant velocity (2p — 1)(1 — p + — p_), as is easily 
determined from conservation of mass. It is then natural to ask what this discontinuity of 
the field u means at the microscopic scale. 

At the microscopic level, a configuration of the ASEP at time t is fully specified by the 
occupation numbers Tj(t), taking values and 1, of all the sites on the lattice. To describe 
the profile of the shock, we must first locate the position of the shock in each configuration. 
This is a nontrivial problem, since we have to distinguish between the variations of density 
due to intrinsic fluctuations and those due to the presence of the shock. Fortunately there 
is, for this system, a simple way of defining the position of the shock at the microscopic 
level: the introduction of a single second class particle into the system. (Alternative and 



more general methods will be discusssed in [13].) The second class particle acts like a hole 
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1 2 - 


-> 2 1 
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with respect to the original, first class particles and like a particle with respect to holes, and 
thus its presence does not affect at all the dynamics of the first class particles. Specifically, 
during any infinitesimal time interval dt, exchanges occur on each bond as follows: 



with probability pdt , 

(1.2) 

with probability qdt , 



where a 0, 1 or 2 on a site means that this site is occupied by a hole, a regular (first class) 
particle or the second class particle. 

It can be shown jOJ] (see also [fl6|, [I7|, and [l|] for a heuristic argument) that the velocity 
of the second class particle in a uniform environment of first class particles at density p is 
(p ~ <?)(! — 2p), so that far to the left of the shock the second class particle moves at a 
velocity (p — q)(l — 2p_), faster than the shock velocity (p — q)(l — p- — p + ), whereas far 
to the right of the shock the second class particle has a velocity (p — q)(l — 2p + ), slower 
than that of the shock. Consequently, the second class particle is attracted to the shock 
and can serve as a marker for its position. It has in fact been proved |15| that the second 



class particle moves with velocity (p — g)(l — p_ — p + ) and |TJ, [L5| (see also |14|], fl9fl-f22| 
that there is an invariant measure for the system viewed from this second class particle, in 
which the asymptotic densities are p±. 
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In the present paper we describe exactly this invariant measure. Our approach is an 
extension of a matrix method which has been used in a number of situations [p4|H35|L in 
which the weight of each configuration is written as the matrix element of a matrix product. 
As in 0, there are three possibilities for each matrix in the product, D, A and E, depending 
on whether the corresponding site is occupied by a first class particle, occupied by a second 
class particle, or is empty. We show in Section 2 and Appendix A that when the matrices 
D, A and E satisfy certain algebraic rules, these matrix products furnish weights for the 
invariant measure we seek. In contrast to most previous cases in which this method was 
used, however, here we write probabilities of events in the invariant measure, using matrix 
products, directly in the infinite system. This approach was used in [3I to recover the 
results of |IJ. 

The invariant measures are parametrized by the pair of densities p_ and p + . There is 
a special value of the asymmetry parameter p, given by p/ (1 —p) = — p_)/p_(l — p+), 
at which the measures are Bernoulli, with density p_ to the left of the second class particle 
and p + to its right. This is described in Section 3, where we also show that the algebraic 
rules given in Section 2 imply certain symmetry properties of the invariant measure. 

The full expression for the shock profile, as seen from the second class particle, is 
described in Section 4, and is derived in Section 5 by constructing an explicit representation 
of the algebra of Section 2. 

We study in Section 6 the asymptotic behavior of this profile at large distances from 
the second class particle, and show that it decays at an exponential rate to p_ or p + . The 
decay length is a function of the parameters p and p± which, however, becomes independent 



of the asymmetry p for p/(l — p) > y p+(l — p~)/ p~(l — p+). 

Finally, in Section 7 we show that in the weak asymmetry regime, where p = (1 + e)/2 
with < e -C 1, the shock profile as seen from the second class particle can be understood as 
the convolution of the hyperbolic tangent profile predicted by the viscous Burgers equation, 

U t + [U(l - U)}y = (l/2)Uyy, (1.3) 

which is known [J7J to describe this regime (on a longer time scale than that on which 
( |1.1|) holds), with the density of the position of the second class particle in this tanh profile, 
given by the derivative of the profile. 

Throughout the paper, the shock will be characterized by the two asymptotic densities 
p_ to the left of the second class particle and p + to the right of the second class particle, 
which satisfy < p_ < p+ < 1, and by the hopping rates p to the right and q = 1 — p to 
the left, as in ( |1.2|) . We will often express our results in terms of p + ,p_,p and q through 
the parameters 

q 

x=-, a = p + (l-p_), b = p-{l-p+). (1.4) 
p 
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2. The algebra for the partially asymmetric shock measure 



In analogy with |l|, [24] we write the probability of a configuration specified by the 
occupation numbers Tj (tj = 0, 1) of m consecutive sites to the left of the second class 
particle (which is located at the origin) and n consecutive sites to its right as a matrix 
element of the form 



-i 



M J] ^ D + C 1 - T i) E \ A II M + C 1 - T J)E\ i«>- (2- 1 ) 

{i=-m ) {j=l ) 

Here, as in |I], a first class particle is represented by a matrix D, a hole by a matrix E, 
and the second class particle by a matrix A. For example, the probability of finding the 
configuration 1 1 to the left and 1 1 to the right of the second class particle, that 
is, of the configuration 

1 1 2 1 1 0, (2.2) 

is given by 

(w\ D E D A E D 2 E 2 \v). (2.3) 

In Appendix A we show, by an extension of the proof which was given in [0, that if 
the matrices D, E and A and the vectors (w\ and \v) satisfy certain algebraic conditions 



then the weights ( |2.1|) are nonnegative and define an invariant measure for the ASEP, as 
seen from a single second class particle, with dynamics specified by ( |1.2|) . The algebra is 



q)[(l-p-)(l-p + )D + p-p + E\, (2.4) 
q)(l-p_)(l- p+ )A, (2.5) 
q) P+P -A, (2.6) 

(2.7) 
(2.8) 
(2.9) 

At this stage the numbers p + and p_ which appear in (|2.4j )- (|2.9| ) are arbitrary parameters 
satisfying < p_ < p+ < 1. However, we will show in Section 6 that they are in fact 
the two asymptotic densities which are reached as one moves away from the second class 
particle. 

As a consequence of ( 2.1|) and of (2.4)— (|2T9|) the microscopic profile, defined here as 



pDE-qED = 


(P- 


pAE - qEA = 


(P- 


pDA - qAD = 


(P- 


(D + E)\v) = 


l*>, 


(w\(D + E) = 




(w\A\v) = 


I. 



the average occupation (r n ) at position n, is given to the right of the second class particle 
by 

(r n ) = (w\ A (D + E) n ~ l D \v), n > 1, (2.10) 
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and to the left of the second class particle by 



( r _ n ) = ( W \D (D + E) n - L A\v), n < -1. (2.11) 

In contrast to the approach to the shock problem taken in [|IJ], expressions like ( |2.10| ) and 
( [2. 11| ) are valid directly for the infinite system; thus here we avoid completely the difficulty 



of taking the infinite volume limit. 



By an argument similar to that of Sandow |25| we may verify that ( |2.4] )-( |2"1]| ) suffice 
to determine any matrix element containing precisely one factor of A. To do so we will 
show that any matrix element of a product of n + 1 operators can be reduced to a sum of 
matrix elements of products of n operators. Consider, for example, a matrix element of the 
form (w\DO n \v), where O n is a product of n operators which contains a single A. Then 

(w\DO n \v) = x k+1 (w\O n D\v) + l.o.t. 

= -x k+1 (w\O n E\v) + l.o.t. 

= -x n+1 {w\EO n \v) +l.o.t. 

= x n+1 (w\DO n \v) + \.o.t. . (2.12) 

Here l.o.t. (lower order terms) denotes matrix elements of products of n matrices, x = q/p 
as in ( |1.4| ), and k is the number of factors of E in O n . To obtain ( |2.12j ) we have used first 
Q and ([[§, then ([£70, then O and O, and finally (fZ|). Since x < 1, equation 
( |2.12| ) can be solved for (w\DO n \v). This reduction permits the calculation of matrix 



elements of arbitrary length. 

It is important to realize, however, that the algebraic rules do not allow us to calculate 
matrix elements of products of operators which contain either no operator A or more than 
one such operator. For example, (w\DE\v) or even (w \v) cannot be calculated from these 
rules alone. Of course, if one has a representation of the matrices D, A, and E and of 
the vectors (w\ and \v) which satisfies ( [2.4H2.9| ) then one can in principle calculate these 
other matrix elements. The values thus obtained depend on the representation and do not 
appear relevant for the problem we consider here. 

When n is small the reduction described above is easy to carry out and one can 
calculate directly the average occupation numbers of the sites closest to the second class 
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particle, as well as other simple correlation functions. Thus for example 
P+ + P- ~ (1 - %)P+P- 



(n) 



1 + x 



(r-i> = fT ^ , (2.13) 



(nr 2 ) 
(r-m) 

( r -2^-l) 



P+ + P+P- + P~ (-1 - a; + 2x 2 )p + p_(p + + p_) + x(l - x) 2 (p + p_) 2 

1 + x + x 2 (1 + 2 + a; 2 ) (1 + x) 

x(p 2 + + P+P- + P-) + (1 ~ ^) 2 P+P-(P+ + P- ~ P+P-) 

1 + X + X 2 

x 2 (p 2 + + P+P- + P_) (2s - x 2 - x 3 )p+p_(p + + g_) + (1 - x) 2 (p + p_) 2 
1 + X + X 2 (1+X + x 2 ) (1 + x) 



In principle one could compute arbitrary correlation functions in this way, but the 
calculation quickly becomes impractical as the number of sites involved increases. We 
note that, in contrast with the totally asymmetric case x = 0, discussed in there are 
in general correlations between the left and the right of the second class particle. The 
expressions in ( |2.13 ) are symmetric polynomial functions of p + and p_. The symmetry 
follows from (|2.4| )- (|2.6|) , but we could not derive it by an elementary argument from the 
dynamics. 

Finally, we observe that if one has a representation of D, E, \v), and (w\ which satisfies 
( p.4| ), ( |2.7| ), and ( |2.8| ) and for which (w\DE — ED\v) is finite and nonzero, then one may 
satisfy (|2J)-(^) by defining 

A = c (DE — ED), (2.14) 



where c is a constant fixed by the condition ( |2.9|) . The representation constructed 
in Section 5 is obtained in this way. We must emphasize, however, that there are 
representations of interest in which ( [2.14j ) is not satisfied. One such representation is 
used in Section 3 below to study the system at a special value x* of x; note also that if one 
wished to use the algebra ( |2.4} )-( |2T9| ) when x = 1 then ( |2.4| ) would imply that D and E 
commute, so that ( |2.14| ) would be inconsistent with the normalization (|2.9| ). More generally, 
( [2.7D and ( |2.8| ) imply that (w\DE\v) = (w\ED\v), so that (w\DE — ED\v ) can be nonzero 
only if (w\DE\v) and (w\ED\v) are infinite. Thus finite dimensional representations of the 
algebra cannot satisfy ( p.!4[) ; the representation constructed in Section 5 is, as expected 
from this remark, infinite dimensional. 
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3. Elementary consequences of the algebra 



We present here two simple consequences of the representation of the invariant measure 
described in Section 2. 

First we note that there is a special value x* of the ratio x = q/p for which the measure 
becomes a Bernoulli measure with density p + to the right of the second class particle and 
density p_ to the left of the second class particle: 

p+Q-p-) a 

This can be seen by verifying that when x = x* the following formulas define a two- 
dimensional representation the algebra (|2.4|)-(p^9|): 



H = (o,i), I«> = (J) • ( 3 - 2 ) 

It is also easy to check directly that this Bernoulli measure is stationary. When x = x* the 
profile on both sides of the second class particle is flat ((r n ) = p + and (r_ n ) = p_ for all 
n > 1) and the occupation numbers at all sites are independent: 

(r_ m . . . r^n . . . r n ) = (r_ m ) . . . (r^fa) . . . (r n ) = p™p™. (3.3) 

The value ( |3.1|) of x* plays the role of a disorder line in equilibrium statistical mechanics 
models |^ gCj. 

Second, we obtain a number of identities satisfied by the invariant measure for all 
x, p_, and p + . These will be derived under the assumption that A has the special form 
Q2.14D , but this is not a restriction on their validity, since we know that there exists at least 



one representation, constructed in Section 5, for which this is true. The simplest of these 
identities is 

(r„)-(r n+1 ) = (w\A(D + E) n -\DE-ED)\v) 

= (w\{DE-ED){D + E) n ~ 1 A\v) (3.4) 
= (r-n-x) - (r-n) , 

which implies that the shock profile has the symmetry 

(t„) + (r- n ) = (r n+ x) + (T- n -i) . (3.5) 
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From the known results ( |2.13|) for the first sites it follows that the common value is 

(r n ) + (r- n ) = P+ + P- ■ (3.6) 



Note that the right hand side of ( |3.6| ) is consistent with the asymptotic values (r n ) — > p + 
and (r- n ) — > p- as n — > oo. 

One can obtain in the same manner identities involving higher correlations. For 
example, from ( |2.14| ) one has 



(w\A(D + E) n '{DE - ED){D + E) m ' D\v) = (w\(DE - ED)(D + E) n ' A(D + E) m 'D\v), 
(w\A{D + E) n 'D(D + E) m '(DE - ED)\v) = (w\{DE - ED)(D + E) n ' D(D + E) rn 'A\v), 
(w\D(D + E) n 'A(D + E) m '(DE - ED)\v) = (w\D{D + E) n ' (DE - ED)(D + E) m 'A\v), 



(3.7) 



and by choosing n' = n — 1 and m! = m — 1 one obtains symmetry relations involving pair 
correlation functions 

\Vn+m+l) v rH-lTn+rri+l) V — n—lTm) Y—n^m/i 

{^~n^"n+m) (Vn+m+l) V — n— m— rn) V— n— mT— rn) > (3-8) 

\J~—nT~rri) v — n^m+l) v — n— m— 1^— m— l) v— n— m— l^*— m) j 

which, through linear combinations, lead to the fact that (r n T n+m ) + (r_ n r m ) + (r_ n _ m r_ m ) 
does not depend on n or m and therefore is given from ( |2.13| ) by 

> =pl + P+ p_ + p 2 _ (3.9) 

Note that this common value is again consistent with the asymptotic densities p± and the 
independence of the occupation numbers r± n , r± m as n,m — > oo (see Remark 6.1). 
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4. The expression for the density profile 



Here we present an expression for the profile (r n ), based on an explicit representation 
of the algebra ( p.4| )-( pl9|) which will be constructed in Section 5. The formula for the profile 
to the right of the second class particle is 

oo 

(r n )=p + + £ Pn-i(k)F k , forn>l, (4.1) 

k=~ oo 

where the P n {k) are defined by the recursion 

Po{k) = S kfi , (4.2) 
P n+1 {k) = aP n {k-l) + {l-a-b)P n {k)+bP n {k + l), (4.3) 

and F k is given by 

F k = —^(a 2 f k+2 ~ a{a + b)f k+1 + b(a + b)f k ^ - b 2 f k _ 2 ), (4.4) 



with 



x 



fk = k r -^, for k ± 0. (4.5) 



fo can be defined arbitrarily because its contribution to ( [4.1|) cancels out: indeed, it is 
easy to check that ( |4.2|) and (|4.3| ) imply that a k P n (—k) = b h P n (k) for all k, so that the 
coefficient of fo in ( |4.1| ) is always zero. (We will introduce below a particular value of fo 
which is convenient in intermediate computations.) The profile to the left of the second 
class particle — (r n ) for n < — 1 — can easily be obtained from the symmetry ( |3.6| ). 

More complicated correlation functions have similar expressions, which can be derived 
in the same way or, in some cases, directly from ( |4 . 1| ) and the algebra. For example, from 
(^.4|)-(^9|), one can show that 



(l-x)D 2 \v) = [D(D + E)-x(D + E)D-(l-x)(l-p + -pJ)D-(l-x)p+pJ\\v), (4.6) 
and this implies that 

(1 - x)(r n r n+ i) = (r n ) - x(r n+ i) - (1 - x)(l - p+ - p-)(r n ) - (1 - x)p+p-, (4.7) 
which gives, using ( [4. 1|) , 

oo 

(r n r n+ i) = p 2 + + Pn-i{k)G k , (4.8) 

k=— oo 
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where 

G k = (p+ + p-)F h + ^—[{a + b)F k - aF k+l - bF k ^] . (4.9) 

1 — x 

Remark 4.1: The quantity P n (k) given by (|4.2| ) and (|4.3|) can be interpreted as the 
probability of finding a biased random walker on site k at time n, given that it was at the 
origin at time 0. We did not find a simple physical interpretation for the presence of this 
random walk in the expression (|4.1|) for the profile but we believe that understanding the 
origin of this biased random walk would give a better insight in the whole problem of the 
description of a shock as seen from a second class particle. 



5. Representation of the algebra 

We now describe the explicit representation of the algebra (|2.4j)-(]2"!]|) used to derive the 



formulas of Section 4 for the profile and to study the weak asymmetric limit in Section 7. 
The operators D, E, and A will act on an infinite dimensional vector space with basis 
{\n) | n = 0, ±1, ±2, . . . }. Let us define the left shift operator L, the two diagonal 
operators 5* and T, and a vector \v) and dual vector (w\ by 

L\n) = \n-l), (5.1) 
T\n) = x n (l +x n )- 1 \n) , (5.2) 
S\n) = [x n (l + X*)- 1 - x^il + x"- 1 )- 1 ]^) , (5.3) 

oo oo 

\v) = £ |n>, (w\ = H • ( 5 - 4 ) 



n=— oo 



These are easily seen to satisfy the following relations: 

pLT - qTL = (p - q)TLT, V TL~ X - qL~ l T = (p - q)TL~ x T\ (5.5) 

(p - q)T(aL - 6L _1 )T = paLT - qaTL - pbTL' 1 + qbL^T; (5.6) 

TL-LT = -LS, L- l T - TL' 1 = -SL" 1 ; (5.7) 

L\v) = \v), (w\L=(w\. (5.8) 

Next we define the operators of our representation: 

D = p^p + + aL-(y/^L-Vb)T(y/^+VbL- 1 ), (5.9) 

E = (l-p_)(l-p+) + 6L- 1 + ( v ^L- v / 6)T( v ^+v / 6L- 1 ), (5.10) 

A = -(a-b)- 2 (DE - ED). (5.11) 
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Then the formula (|2.4j) for pDE — qED follows immediately from (|5.6|) , and as noted in 
Section 2 the relations ( p.5| )-( ^6| ) for pAE — qEA and pDA — qAD are automatically 
satisfied by (|5.11|) . Moreover 



D + E = aL + (1 -a-b)I + bL~ 1 , (5.12) 
so that (D + E)\v) = \v) and (w\(D + E) = (w\. Finally, ( [5.7| ) implies that 

A = -(a - b)- 2 {^L - Vb){aLS - bSL- 1 )^ + VbL- 1 ), (5.13) 

and from (w\S\v) = —1 it follows that = 1. 

Note that, as expected from the remark following ( [2.14|) , the matrix elements {w\v), 
(w\DE\v) and (w\ED\v) are all infinite. In this representation the presence of a factor A 
in a matrix product — or more specifically, that of S (see (|5.13| )), a diagonal matrix whose 
diagonal elements Sjj decrease exponentially fast to as j ' ^ ±oo — renders the (w\ ■ \v) 
matrix element of that product finite. 

We now turn to the evaluation of (r n ). First we observe that, for k > 0, 



\ x j-k 



£ T 

j=-N X 



N 



x j-k I _j_ x j 



1 - X k [ 1 



X J 



x- 



j-k 



j=-N 
N 



£ T 



X J 



+ 



X 



1 + xi~ k 

N 



. £..T 



, _|_ x-? 1 — T k \ ^ — 4 1 -I- t3 

j=-N x T x 1 x W=JV-Jfe+l ^ x j=-N-k 



-N-l 



(5.14) 



From ( 5.14 ) we have, for fc > 0, 



HSL fc 2» = lim V - . . , , . . , 



lim 

TV— »oo 



A' 



1 £ fe \ j =J 0Cfc + i 1 + x3 j=tjtf-k 1 + xJ 



fe+i 



1 - a; fc+1 + 

\3=N—k 



N j -N-l 

y - e - 

V , 1 + z J . tr" , , 1 



j=-N-k-l 



X 3 



— fk+1 ~ fki 

with fk = kx k /(l — x k ) as in (4.5). From the easily verified identity 

(iv\SL j T\v) + {w\SL~i- l T\v) = -1, 



(5.15) 



(5.16) 
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it then follows that ( |5.15|) also holds for k < — 1, with f k again given by ( |4.5| ). Now let 
us define fo = fi — (w\ST\v); then (|5.15|) holds for k = and also k — —I (using Q5.16P 
again), and thus for all k. It is in fact not necessary to evaluate fo, by the remark following 

(El- 

Using (|5.9j), ( [5. 13| ) , and ( |5.15| ) we have, for all k, 

(w\AL k D\v) = p^p + + a+{a-b)- 1 (w\(aLS -bSL^)L k {aL -bL~^)T\v) 
= p+ + {a- by l (w\S(a 2 L - ab - abL' 1 + b 2 L- 2 )L k T\v) 
= p+ + (a- b)- l (a 2 f k+2 - a{a + b)f k+1 + b{a + b)f k ^ - b 2 f k „ 2 ). (5.17) 

The formula ( |5.12 ) for D + E yields by induction 

(D + E) n = ^2 Pn{k)L k (5.18) 

k 

for n > 0, where the P n {k) are defined by (|4.2| ) and ( |4.3|) . Thus for n > 1, 
(r n ) = ^^(.D + E)"- 1 /}^) 

oo 

= ]T P^^A;)^!^^^) (5.19) 

fc=— oo 

oo 

= P+ + ^ZTl P n-i{k){a 2 f k+ 2-a{a + b)f k+l + b(a + b)f k _ l -b 2 f k ^ 2 ), 

and this completes the derivation of the formula ( [4.1|) for (r n ). 

From the expression (|4.1|) it is possible to evaluate various asymptotics of the profile. 
The two main limits that one might wish to consider are the limit n — > oo describing the 
tail of the profile, discussed in the following section (where we also comment briefly on the 
p + \ p_ limit), and the limit of a weak asymmetry, x — > 1, discussed in Section 7. 



6. Spatial asymptotics of the density profile 

We now discuss the behavior of the profile (r n ) at large n, basing our discussion on 
( f4.1|) . Because a > b, the biased random walker governed by (|4.3|) goes to +oo asn-> oo, 
with distribution P n (k) concentrated around k ~ n(a — b). Moreover, F k approaches as 
k approaches oo, decaying exponentially (as kx k ), so that the sum in ( [4.1| ) vanishes in the 
n —>■ oo limit and thus the asymptotic density on the right is p + , i.e., 

( T n) — P+ as n — >• oo. (6.1) 
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From the symmetry ( |3.6| ), the asymptotic density at the left is p_. 

To obtain the approach of (r re ) to its n — > oo limit, we need to estimate how the sum 
on the right hand side of (O) vanishes. As observed above, decays exponentially as 
k — > oo, and it is easy to check that it approaches the constant value —(a — b) as k — > — oo. 
Because of this behavior there are two cases to consider: the sum is dominated either by 
large values of k or by values of k close to 0, depending on the value of x. 

Case 1: x> (b/a) 1 / 2 

We first explore the consequences of assuming that the sum in (|4.1|) is dominated by 
large values of k. If this is so, one can use the approximation 



and then the identity 



k 



fk - kx\ (6.2) 



oo / h \ n 

x k P n (k) = fl -a-b + ax + -J . (6.3) 



This leads to 



(T n )- P+ ^(n-l) iaX - b/x)2{aX - a - b + b/x) (l-a-b + ax + h X 2 

a — o \ x J 

2a 2 x 2 - (a + b)ax - (a + b)b/x + 2b 2 /x 2 ( , b \ . . 

+ y -—J \ J ' / —[l-a-b + ax+-\ . 6.4) 

a — b \ xl 



The values of k which dominate the sum (^3) are k ~ n(ax — b/x)/ (1 — a — b + ax + b/x), so 
that the assumption that large values of k dominate the sum in (|4.1|) would be inconsistent 
if x 2 < b/a. Conversely, when x 2 > b/a, (|6.4j) gives the asymptotic behavior of (r n ). 

Case II: x < {b/a) 1 / 2 

Since is constant for k — > — oo and P n (k) increases with k for k < 0, it is clear that 
the sum over k must be dominated by values of k near k = (i.e. which do not scale with 
n). It is convenient to use the following expression for the solution P n (k) of (|4.3|): 



Pn(k) =£"(!-«- Z ft^a^. (6.5) 



For large m and for of order 1, one can write 



m 

m+k 
v 2 / 

m 

m+fc 



2 m \/ — ~ (l — — ) , for m + even, 
V nm \ 2m J 

0, for m + k odd. (6.6) 
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Then 



l/2^/7i-lV n „ mWnW2 2 m ~ fa\ k / 2 ( A; 2 \ 



1 2 ^ fn-V 



E ('^1( 1 - a - & r m " 2m (^) m/2 TT-^- (6.7) 



2m 3/2 

Here we have used the fact that Y,'k=-cc,{ a /b) k/2 Fk = (see (fPl ) and flUf)) and set 



fc=-oo VU/ U U fc=-oo 1 X KU/ 

The extra factor of 1/2 in (|6.7| ) comes from the fact that we have replaced a sum in which 
m and k have the same parity by a sum over all m and fc. Finally, using the fact that 
(^x n ~ m y m m~ a ~ (x + y) n+a (ny)~ a , since this sum is dominated by values of m near 
ny/(x + y), (|6.7|) becomes 



By pairing the ±fc terms in (|6.8| ) we see that if x < x* = b/a then Z > and 
asymptotically (r n ) > p + , while if x > x*, asymptotically (r n ) < p + . Numerical 
computation indicates that (r n ) — p + has the same sign for all n and is in fact monotonically 
decreasing for x < x* and monotonically increasing for x > x*. Note also that Z vanishes 
for expected from the special nature of the measure at this value of x (see 

Section 3). 

When x = 0, (|6~8|) becomes 



Aab 



(6.10) 



leading to 

^-^- 2^-^ ^ ll -^ a - Vt>rH - (6 - U) 
This agrees with (6.26) of H (and with (7.6) of that paper after the correction of a misprint: 
the exponent 2i — 1 should be 2i + 1). 

By a modification of the above calculation one may also find the asymptotics when 



p + = P-. The result is independent of x and is that given in [|TJ: (r n ) ~ p + y p(l — p)/nn 
when n>l, where p = p + = p_. 



14 



The various asymptotics for the profile derived here are summarized in Figure 1 . Figure 
2 shows the profile for some typical parameter values. 

Remark 6.1 Using the explicit representation of the algebra developed in Section 5 it can 
be verified that the measure is asymptotically Bernoulli as n —>■ ±oo, i.e., that if T is any 
operator product containing j factors of D and k factors of E then 

lim (w\A(D + E) n T\v) = p> + (\ - P+ ) k , (6.12) 

with a similar limit at — oo. The approach is exponentially fast, with the same decay rate 
as for the profile, that is, (1 — a — b + ax + b/x) n when x > (b/a) 1 ^ 2 and [1 — (yfa — Vb) 2 )] n 
when x < (5/a) 1//2 . More generally, if A is a product of D's, £"s and precisely one A, then 
(w\ A(D + E) n T\v) — > A|t>)p+(1 — p + ) k asn^ oo, with the same rate of approach (and 
again a similar limit at — oo). 



7. The weakly asymmetric limit 



In this section we study the limiting behavior of the model as the asymmetry e = p — q 
becomes small, with p± fixed. In this case the profile (r n ) becomes very broad and as e — > 
depends only on the scaled variable y = ne. To verify this, note that for e small and k large 
it follows from 

x = q /p= l _2e + 0(e 2 ) 

that 



f k+p = (k + p)- 



x 



k+p 



X 



so that 



(a -b) 



1 — x } 

k 



+ p 



1 — X 1 



— 2pke- 



x k ) 2 



(7.1) 
(7.2) 



x 



2ke- 



(7.3) 



1 — x k (1 — X" 

Now for large n the probability P n {k) is concentrated around k = n(a — b); if we introduce 
the variable y = ne, use the approximation (|7.3|) , and set k = (a — b)y/e we obtain 



(r„> 



P+ + P- 



1 - e-^v 



-4\y 



+ A 



coth2Ay 



1 _ e -iXy)2 
sinh 2 2 Xy 



(7.4) 



where A = (a — b)/2 = (p + — P-)/2. It is easy to see that the terms neglected in ( |7.4j ) 
are 0(e) when y is of order 1. Note that the asymptotic behavior of (|7.4|) agrees with the 
x — > 1 limit of 
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This scaling form for the profile can be understood in terms of the viscous Burgers 
equation for a hydrodynamic density variable u(y, s): 



U S + [U(l - U)]y = (l/2)Uyy. (7.5) 



It is shown in |$6 1 , |37| ( see a l so HH ) th & t the weakly asymmetric macroscopic limit of the 
ASEP, in scaled variables y = ne as above and s, related to the microscopic time t by s = e 2 t, 
is indeed (|7.5| ). For example, if the initial state of the particle system with parameter e 
is a product measure with density u (en) at site n, then the observed density u(y,s) at 
position e~ l y and time e~ 2 s is the solution of ([7.5|) with initial condition u(y, 0) = uo(y). 
Equation (|7.5|) is diffusive and does not exhibit discontinuous shocks; instead of the jump 
in density from p_ to p + we now have a smooth profile u(y, s) = p(y — vs), traveling with 
velocity v — 1 — p + — p_, with 

P(y) = (P- + P+)/2 + Atanh2Ay, (7.6) 

and A = (p + — p~)/2. Thus the shock width on the microscopic scale diverges in this limit 
as e _1 . 

We can now interpret (|7.4j) , after putting y = en, as the shock profile seen from 
a second class particle which does not have a fixed location relative to the tanh profile; 
instead, its position is distributed in such a way that it sees all densities between p_ and 
p + with equal probability. Then the probability of finding the second class particle between 
y and y + dy is dp(y) = Q(y) dy, where 

Q(y) = p'{y) = A cosh" 2 2Ay, (7.7) 

P+ ~ P- 

and the limit of the one particle density function will be 

2Xy 



f°° P+ + P- 
lim(r Le -i 2/J ) = / p(y + z)Q(z) dz = + A 

e J —CO Z 



cothAy —n — — 

sinh 2Xy 



(7- 



in agreement with (|7.4j ); here [^J denotes the greatest integer k satisfying k < z. 

Remark 7.1: In [IIS] the shock is viewed, not from the position of a second class particle, 
but from a location defined in a rather less direct manner. With that definition the profile 
seen in the weakly asymmetric limit is p(y), not the convolution ( |7.8| ). 

There are various possible direct justifications for the density (|7.7|) . A heuristic 
argument may be based on a comparison of two systems: in one of these the second class 
particle is treated as a first class particle and in the other as a hole. In the latter case 



16 



the shock position, however determined, is effectively translated a macroscopic distance 
e/(p+ — p~) relative to its position in the former. Thus the density of the second class 
particle is (p(y + e/(p+ — — p(y))/e ~ p'(y)/(p + — pJ). A rigorous version of this 
argument is given in [[15]. 

The above picture is confirmed if we compute other aspects of the e — > limiting 
behavior of the entire invariant measure. Consider first the distribution of occupation 
numbers at a finite number of sites with both microscopic and macroscopic spacing. To 
give a concrete example we consider 

(w\A(D + E) ni DED(D + E) n2 ~ ni EEDD\v) (7.9) 

where n\ = [j/i/ e J an d n 2 = L2/2AJ f° r some macroscopic positions y±, y 2 with < y\ < y 2 ] 
this is the probability that of three specific sites at y\ the first and third are occupied and 
the second empty, and that the first two of four sites at y 2 are empty and the second two 
occupied. Then 

lim (w\A(D + E) ni DED(D + E) n ' 2 ~ ni EEDD\v) 

/oo 
p{y x + z)\\ - p( Vl + z))p(y 2 + z)\l - p(y 2 + z)f Q(z) dz. (7.10) 
-oo 

We may think of the right hand side of ( |7.10|) as arising from a distribution of particles in 
which sites at macroscopic positions y are occupied with probability p(y) and all occupation 
numbers are independent; this distribution is viewed from a random position z (the position 
of the second class particle) distributed according to Q(z) dz. Equation (|7.10|) — in fact, its 
generalization to an arbitrary number of sites — will be proved in Appendix B; in particular, 
this furnishes an alternate derivation of the special case ( |7.4| ) of one site. 

We can also show that, in a certain sense, the typical configuration in the e — ► limit 
has on the macroscopic scale the shape of the tanh profile (|7.6|) . To do so we introduce 
a coarse graining on an intermediate scale e _/3 , where < /3 < 1: for each macroscopic 
position y we define a random variable B e (y), the empirical density in a block of size e _/3 
at y, by 

B e (y) = ^ E nv/el+k- (7.H) 
fc=i 

Now consider, for very small e, a configuration of the system which after coarse graining 
looks like the tanh profile ( |7.6|) on the macroscopic scale. The second class particle will 
be located at some macroscopic position z relative to the center of this tanh profile and 
hence will see locally a density p(z); thus B e (0) ~ p(z) and B e (0) determines the position 



17 



z by z ~ p _1 (i? e (0)). At position y relative to the second class particle the density will be 
p(z + y), and we are led to expect a strong correlation between the empirical densities at 
positions and y relative to the second class particle, given by 



B e {y)~p(y + p -\B e {0))). 



(7.12) 



We show in Appendix B that ( |7.12| ) holds with arbitrary accuracy, and with probability 
arbitrarily close to one, for sufficiently small e, that is, that for any rj > 0, 



lim Prob 



B £ (y)-p(y + p- L (B € (0)))> V 







(7.13) 



(convergence in probability). Equation ( [7.1 3| ) essentially says that, with probability one, 
the coarse grained densities in the e — > limit lie on some translate of the profile p(x). 



8. Concluding remarks 

We have found a family of stationary measures for the ASEP as seen from a second 
class particle, parametrized by two numbers p_ and p + , with < p_ < p + < 1, which 
correspond to densities of the asymptotic Bernoulli measures at ±oo. The results, which 
hold for all p > 1/2, generalize those of |1[, in which such measures were obtained for 
the fully asymmetric case p — 1. Our derivation, in contrast to that of previous works, is 
carried out directly in the infinite system. 

For a certain value of the asymmetry, x = x*, we found a two dimensional 
representation of the algebra (see ( |3.1| ) and ( |3.2| )). In fact, it can be shown that a 
2r dimensional representation exists when x r = x*, r = 1,2,.... Finite dimensional 



representations of other ASEP algebras have been found in |29|, [43]. 

In the weak asymmetry limit p — > 1/2, the measure can be understood in terms of the 
solutions of the macroscopic viscous Burgers equation. Several recent works have studied 
the gap of the evolution operator or generator f4T| , 42], as well as the diffusion constant [32]. 
on a ring of size L, in the double limit p ^ 1/2, L ^ oo. A simple scaling form was found 
for the diffusion constant; it would be nice to see whether a scaling form in either problem 
could also be derived from a macroscopic equation. 

The microscopic shock profile derived in this work contains both intrinsic features of 
the shock and the fluctuations of the position of the second class particle. One can imagine 
many other ways of locating the shock with an alternate definition, the expression 



for the profile would certainly be different. However, one expects that there exist intrinsic 
properties of the shock which are independent of the definition of its location. One class of 
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such properties would be the values of sums of expectations of all translates of quantities 
whose expectation values vanish in the asymptotic Bernoulli measures with densitites p + 
and p_, e.g., 

oo 

£ (( T i ~ p-)Ti + n(r i+ n+m ~ P+))- (8.1) 

i=— oo 

We have evaluated explicitly such sums for the totally asymmetric case p — 1, using the 
invariant measure as seen from the second class particle, and will report on this work in a 
later publication [HJ. From this point of view, identities such as ( |3.4| ) reflect the fact that 
the second class particle can be considered, in the sum, as either a first class particle or a 
hole. 
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Appendix A. Verification that the algebra yields an invariant measure 

In this appendix we show that if the operators D, E, and A and the vectors \v) and 
(w\ satisfy the algebra ( |2 . 4|) — ( p79|) then the formula ( |2.1|) yields an invariant measure for 
the ASEP as seen from a second class particle. We begin by writing down the conditions 
which must be verified. The algebra allows us to calculate directly the probability of 
finding configurations r_ m , . . . , r_i to the left, and t±, . . . , r n to the right, of the second 
class particle, where the r, are thought of as elements of the set {0, 1}. For convenience in 
describing the dynamics and consistency with the notation of (|1.2j) we will write this local 
configuration as r = (r_ m , . . . , r_i, r , Ti, . . . , r n ); in this notation the symbol r always 
takes the value r = 2 and represents the presence of a second class particle at the origin. 

Suppose then that for each such r we have defined a weight W{r). Then these weights 
define an invariant probability measure for the ASEP if they satisfy (i) normalization 
conditions: 

W(r ) = 1, 

W(T_ m , . . . ,T_i,T ,Ti, . . . ,T„) = W(T- m ,...,T- 1 ,T ,Ti,...,T n ,a) (A.l) 

(7=0,1 

= J2 W(a,T- m ,...,T- 1 ,T ,n,...,T n ); 

0=0,1 
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(ii) positivity: 

W(r_ m , . . . , r_i, t 0) ri, . . . , r n ) > ; (A.2) 

and (iii) for each configuration r a condition corresponding to the stationarity of the 
probability of r. To write down these stationarity conditions we will use the following 
notation: if r = (r_ m , . . . , r„) is a system configuration and — m < i < n — 1 then r l,l+1 
denotes the configuration obtained from r by a particle jump across the bond between sites 
i and i + 1; if i, % + 1 ^ this is an interchange of r« and r i+ i, but if either i = 0ori + l = 
it is a shift of the configuration relative to the second class particle: 



-0,1 



T J 



1,0 



r 



if J = -1, 
if .7 = 0, 

if — m — 1 < j < —2 or 1 < j < n — 1; 

if J = 1, 
if J = 0, 

if — m + l<j<— lor2<j<n + l; 
if j = i, 

if j = i + 1, when i + 1 ^ 0. 

if j 7^ M + !> 



(A.3) 



Now during some time interval <it there is, for each pair of adjacent sites i, i + 1 with 
7? 7^ T i+ii some probability of exit from the configuration r by an interchange r — » r l,l+1 , 
and some probability of entrance into r by an interchange t 1 ' 1+1 — > r: if Tj Tj + i has the 
form 10, 12, or 2 then the exit probability is pdt and the entrance probability qdt, 
while if Tj r i+ i has the form 1, 2 1, or 2 then these probabilities are reversed (see ( P-.2| )). 
Moreover, there are probabilities for exit from and entrance to r due to exchanges across 
the ends of r, between sites — m and — m — 1 or n and n + 1; if n = or m — 0, these 
exchanges involve the second class particle. Combining these effects we find the equation 



n-l 



n-1 



= J2 ' l-P w ( T ) + qW(r l > l+1 )} + " [~qW(t) + P W(t 1 > 1+1 ] 

i=—m i=—m 

+Xi(r- m )[-qW(0T)+pW((0r)- m - 1 ^ m )} 
+Xo(r^ m )[-pW(lr) + ^((lr)-™- 1 '-™)] 
+Xi(r n )[-pW(r0) + qW((r0) n ' n+1 )] 
+Xo(r n )[-qW{rl) + P W \(rip n+1 )] , 



(A.4) 



where the singly (respectively doubly) primed sum is over indices i for which TjTj+i is 10, 
12 or 20 (respectively 01, 21 or 02), and the indicator functions Xi an d Xo are defined by 



if a 
if a 



1 or a 
0, 



Xo(o-) 



if a 
if (X 



or a 
1. 



(A.5) 
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Equation ( |A.4| ) demands special attention if all the other than tq are the same, 
say Ti = 1, % ^ 0. Then exchanges involving the second class particle will not always 
change the local configuration. For example, if r = (1, 1, 1, 2, 1, 1) then after the exchange 
t — > r ' 1 = (1,1,1,1,2,1) the local configuration is still r if there was originally a 
particle on site 3 (the third site to the right of the origin). Thus ( |A.4j ), which includes 
a term — qW(r) from this exchange, should be corrected by a term +qW(rl). But 
( |A.4|) also contains a term +qW(r~ 1,0 ) from the exchange r -1 ' = (1,1,2,1,1,1) — > r; 



in this case the local configuration was already r if there was before the exchange a 
particle at site —3. Thus ( |A.4| ) should also be corrected by a term — qW ((lr)^ 1 ' ) . Since 
rl = (1t)~ 1,0 , these two corrections cancel. Said otherwise, (|A.4|) counts the transition 



(1, 1, 1, 2, 1, 1, 1) — > (1, 1, 1, 1, 2, 1, 1), which does not change r, twice, with opposite signs. 
The argument is easily seen to be quite general and to apply also to leftward jumps of the 
second class particle and to the case Tj = 0, i 0. We conclude that ( |A.4j ) is correct in all 
cases. 

We will now show that the weights defined by (|2.1|) , 

-1 n 

W (T_ m ,...,T_ 1 ,T ,T 1 ,...,T n ) = (w\ J] [r i D + (l-r l )E}Al[[T J D + (l-T j )E}\v), (A.6) 

i=—m j=l 



satisfy ( |A.1|) , (|A.2|) , and ( |A.4|) and hence provide an invariant measure for the ASEP. ( |A.1| ) 



follows immedately from ( |2.71 )-( |2T9| ). We will verify ( |A.2| ) by induction onm + n, the total 
number of D and E operators in the product. Since (w|A|t>) = 1 by fl2.9p the case m+n = 
is immediate. For general m + n it suffices to show that 

(w\E m AD n \v) > 0, (A.7) 



since by repeated use of the relations (|2.4| )- (|2.6| ) we may express the matrix element of any 
product with m operators E and n operators D as x k (w\E m AD n \v), for some k, plus a 
linear combination, with positive coefficients, of matrix elements of products with a smaller 
number of D or E operators. 

To verify ( | A . 7| ) we define new operators D and E by D = p+p- + D and 
E = (1 — p + )(l — p_) + _E and show, again by induction, that 

(w\E m AD n \v) > 0, (A.8) 



from which ( A.7|) (with strict inequality) follows immediately. The operators D and E 
satisfy 

DE — xED = (1 — x)ab, DA — xAD = 0, AE — xEA = 0, 
(D + E)\v) = (a + b)\v), (w\(D + E) = (w\(a + b), 
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and hence also D k E — x k ED k = ab(l — x k )D k 1 and DE k — x k E k D = ab(l — x k )E k 1 for 
any k > 1. Then an argument as in ( [2.12| ) leads to the recursions 

1 



w\E m+1 A\v) 



1 - x m+2 



(a + b)(l-x m+1 )(w\E m A\v) 

-ab{l - x^iwlE^Alv) 



(A.9) 



w\E m AD n+1 \v) = (a + b)(w\E m AD n \v) - x n+l (w\E m+1 AD n \v) 



n+l , 



-im+1 



-ab(l - x n ){w\E m AD n - 1 \v)- 1 (A.10) 

similar recursions hold with m and n interchanged, hence (w\E m AD n \v ) = (w\E n AD m \v ). 
From (|A.9| ) it follows by induction on m that 



(w\E m A\v) 



i m+1 - b 



,-m+l 



1 



a - b 1 - x m+1 ' 
and then from ( |A.10|) and ( |A.11|) , by induction on n, that for n < m, 



(A.n; 



< w 



E m AD n \v) = Y,( ab Y 

k=0 



im +n+l-2fc _ h m+n+l-2k ^ _ ^ n™ =fc+1 (l - X*] 



n 



n—k i 
3=0 \ 



X 



(A.12) 



(the identity (a + b)(a j - V) = (a j+1 - V +1 ) + afo(a i_1 - V' 1 ) is needed in checking ( |ATl| ) 
and ( |A.12| )). The expression ( |A.12|) is clearly positive, so that ( |A.8|) holds. 

Finally, we must verify ( |A.4|) . To do so, we partition the matrix product in ( |A.6| ) 
into blocks of consecutive identical matrices — that is, blocks of D's or E's, together 
with one block consisting of a single A. When the form ( |A.6| ) is substituted into the 
stationarity condition ( |A.4| ), each resulting term arises from a possible exchange at some 
block boundary — between blocks or at the end or beginning of the product — and contains 
a corresponding factor ±(pDE — qED), ±(pAE — qEA), or ±(pDA — qAD). These factors 
may be simplified with the fundamental relations ( |2.4j )-( p^6|) , which we will use in several 
forms: 



pDE - 


qED 


= (p- 


q)(eD + dE) 






= (p- 


q){{d-e)E + e{D + E)) 






= (p- 


q)((e-d)D + d(D + E)), 


pAE - 


qEA 


= (P- 


q)eA, 


pDA — 


qAD 


= (P- 


q)dA. 


and e = (1 - 


"P-)(l 


-P+)- 





(A.13) 
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The net effect of this simplification on each possible block boundary is summarized in 
the following table: 



Boundary Contribution to ( |A.4|) Simplified contribution 



DE 


(w 


(-pDE + qED) 




v) 


(p- 


-q) 


(w 




(-eD - dE) ■ ■ ■ 


v) 


ED 


(w 


■ ■ ■ (+pDE - qED) 




v) 


(p- 


-q) 


(w 




(+eD + dE) ■ ■ ■ 


v) 


DA 


(w 


■ ■ ■ (—pDA + qAD) 




v) 


(p- 


-q) 


(w 




■ (-dA) ■■■ 


v) 


AD 


(w 


■ ■ ■ (+pDA - qAD) 




v) 


f 

(p- 


-q) 


(w 




■ (+dA) ■ ■ ■ 


v) 


AE 


(w 


■ ■ ■ (~pAE + qEA) 




v) 


(p- 


-q) 


(w 




■ {-eA) ■■■ 


v) 


EA 


(w 


■ ■ ■ (+pAE - qEA) 




v) 


(p- 


-q) 


(w 




■ (+eA) ■■■ 


v) 


(w\D 


(w 


(+pDE - qED) 




v) 


(p- 


-q) 


(w 


((e 


-d)D + d) ■■■ 


v) 


(w\E 


(w 


{-pDE + qED) 




v) 


(p- 


-q) 


(w 


((e 


-d)E-e) ■■■ 


v) 


(w\A 


(w 


{—pDA + qAD + pAE - qEA) ■ ■ ■ 


v) 


(p- 


-q) 


(w 


((e 


- d)A) ■ ■ ■ 


v) 


D\v) 


(w 


(-pDE + qED) 


v) 


(p- 


-q) 


(w 




■ ((d-e)D-d) 


v) 


E\v) 


(w 


(+pDE - 


qED) 


v) 


(p- 


-q) 


(w 




■ {{d-e)E + e) 


v) 


A\v) 


(w 


■ ■ ■ (+pDA - qAD - pAE + 


qEA) 


v) 


(p- 


-q) 


(w 




■■■ ((d-e)A) 


v) 



The last column shows that each block of D's gives rise to two terms, one from each 
boundary of the block, in which one of the factors of D is replaced by the constant d; these 
terms have opposite signs and hence cancel. Two similarly cancelling terms arise from 
each block of £"s. Finally, the left and right ends of the product give rise to additional 
terms in which the original amplitude is multiplied by (e — d) and (d — e), respectively; the 
cancellation of these terms completes the verification of (|A.4|) for the weights ( |A.6| ). 



Appendix B. The invariant measure in the weakly asymmetric limit 

In this appendix we verify the picture given in Section 7 of the weakly asymmetric 
limit of the invariant measure: for e = p — q very small, the measure is approximately a 
convolution of a Bernoulli measure having density profile p(y) with the density Q(y) of the 
position of the second class particle in this profile. Here y = en, so that p(y) and Q(y) (see 
( |7.6| ) and (|7.7|)) vary on the macroscopic scale . 

We begin by describing the key steps in the argument. The probability of occupation 
numbers <7i, <7 2 , • • • cr m at specified sites ri\ < n 2 < ■ ■ ■ < n m is given by 

(w\(aiD + (1 - a 1 )E){D + E) 712 ' 711 - 1 (a 2 D + (1 - a 2 )E){D + E) n3 ~ n2 ~ l ■ ■ ■ A ■ ■ ■ 

...(£> + E) n ™- n —-\o m D + (1 - a m )E)\v) (B.l) 
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We want to compute an approximation to ( |B.1|) that is correct in the e — ^ limit. First, each 
factor (D + E) n in the product has matrix elements (D + E)^ = P n (i—j); since P n (k) is the 
probability distribution of a biased random walker and is concentrated near k = n(a — b), 
we approximate (D + E)"^ 1 ~ (D + E)^ — 5i_j i |_ n ( a -&)j • Next, we approximate each single 



factor of D, E, or A in the product by a diagonal matrix, obtained from equations (|5.9Q , 
(|5.10 ), and ( |5.13| ) by making first the approximation L ~ L^ 1 ~ / and then a simple 
numerical approximation: 



D 



i.i 



P+-2A 



x-> 



1 +xi 



5 



'J 



P 



5L 

2A 



6.. 



y ■ 



A; 



1-P+ + 2X 



1 + xi 



-SjjSij 



2A^ V2A 



With these approximations the inner product ( |B.1| ) becomes 



e n 

j=— oo fc=l . 



^(j + n k )' 
2A 



+ (1 



+n k ) 
2A 



But this is just a Riemann sum for the integral 



2A 



(B.2) 

(B.3) 
(B.4) 



Q(z) dz Yl [vkpiVk + z) + (1 - <7 fc )(l - + z)] , 
fc=i 



(B.6) 



where = en&. Since the integrand vanishes exponentially fast at ±oo there is no difficulty 
with convergence of the Riemann sums. Moreover, we will sketch below a proof that the 
sum of the errors made in the approximations leading to ( |B.5| ) — that is, the difference 
between ( |B.1| ) and ( B.5 ) — vanishes as e — » 0, uniformly in the choices of n k and o k . Thus 
we conclude: if the n k are chosen to depend on e in such a way that lim e ^ en k — Vk then 



limProb{r nfc 



a k | k = 1, . . . , m} 

/OO 771 
Q(z) cfe J] fopfe* + «) + (!- ^)(! " + *)] • (B.7) 
"OO 7. 1 



fe=l 



Equation QB.Tp gives the e — > limiting behavior of the profile, equation ( |7.8[) , and 
the similar behavior ( |7.l0|) in the example discussed in Section 7 for the distribution of the 
occupation numbers of several sites. In general, ( p. 7] ) describes the limiting behavior of 
the distribution of occupation numbers at a finite number of sites with both microscopic 
and macroscopic spacing, since we do not assume that the positions y k are distinct. 
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We now discuss the coarse grained random variables B e (y), the empirical densities on 
an e _/3 scale, where < (5 < 1 (see ( |7.11| )). We want to show that for a typical configuration 
in the e — > limit these variables lie on the hyperbolic tangent profile p(y), after a 
configuration-dependent translation. Due to this random translation, however, the B e 
variables fluctuate even in the e — > limit; to obtain a sharp statement, we show that there 
are no fluctuations in a function of several of these variables, g(B) = g(B € (yi), . . . , B e (y K )), 
when g is independent of translations of the y variables (see (|B.10| ) below). One example of 
this technique is ( |7.12| ). To show the absence of fluctuations in g(B) we want to show that 
(\g(B) — (g(B))\) = 0; the idea then is to use (|B.7| ) to compute expectations of functions 
of the B t (y). 

Specifically, suppose that g(£i, . . . is a function defined for < < 1 and that 
yi,...yx are any real numbers. We will show that 

/oo 
g(p( yi + z),..., p(y K + z)) Q(z) dz. (B.8) 
-oo 

We verify ( |B.8|) first for g a monomial, then a polynomial, then continuous, and finally for g 
bounded and continuous on (p_, p+) K , the case needed in applications. If g is a monomial, 
say of degree N, then we may substitute the definition ( |7.11| ) of B e into g and expand, so 
that (g(B e (yi), . . . , B e (y K ))) becomes a sum of ^"^J^ terms e~ Nl3 (Y[f =1 T mj ). There are at 
most (^e~ < - Ar ~ 1 - )/3 such terms in which not all the r mj are distinct, so that their contribution 
can be ignored, and from ( |B.7| ) the remaining terms give precisely the right hand side of 
( p.8|) . Thus ( p.8|) holds if g is a polynomial and hence, by a uniform approximation 
argument, if g is continuous on [0,1]^. If we apply the latter result, for K = 1, to a 
continuous function hs(£,) satisfying < hs(£) < 1, h$(£) = 1 for p_ + 25 < £ < p + — 25, 
and hs(£) = for ^ < p_ + 5 and £ > p + — 5 we have 

/oo 
h 5 (p(y + z))Q(z)dz, (B.9) 
-oo 

so that Prob{p_ + 5 < B e (y) < p + — 5} — > 1 as 5, e — » 0. Using this result, we may 
extend (|B.8|) to any function g which is bounded and is continuous on (p_, p+) K . To do so, 
we restrict g to (p_ + 5, p + — 5) K and then apply ( |B.8|) to a continuous extension of this 
restriction satisfying the same bound as g, obtaining ( |B.8| ) for g with error which vanishes 
as 5 — > 0. 

Now consider a function ^ as above, bounded and continuous on (p_, p+) K , such that 
for some yi, . . . , yx, g(p(yi), ■ ■ ■ , p{vk)) is translation invariant: for any z, 

g{p{yi), • • • , p(2/*0) = ^(p(?/i + z),..., p(y K + z)). (B.10) 



25 



Then 



lim , ( \g(B e (yi), B e (y K )) - g(p(yi), p{vk)) |) 

e— »U 



g(p(yi + z),..., p(y K + z)) - g(p(yi), . . . , p{y K ))\ Q{z) dz 
= 0, (B.ll) 

so that \im e ^Q g(B e (y l ), . . .,B e (y K )) = g{p{yi), . . . , p(yx)) in probability. 

From (p.ll| ) we obtain the result described at the end of Section 7: if for any y we 
define g{£i,&) = 6 - p(y + and then apply (|B.11|) to g(B e (0), B e (y)), we obtain 

(|7.13|) . Note that g is bounded but is not continuous at £i = £2 = P- or £1 = £ 2 = p + . 



We finally sketch the argument that the errors made in the approximations leading to 
( p.5|) are uniformly small. Control of errors introduced by the approximation (D + E)^ ~ 
bi-j,\n{a-b)\ is straightforward, and details will be omitted. To discuss (p.2| )-( ]B~4|) , we 
simplify the notation by supposing that n\ > . Let us write for the diagonal 
approximation to D introduced in (|B.2| ), and Z)M = uD + (1 — u)D^ for < u < 1, with 
similar notation for E and A. We must then estimate 

f 1 ^{w\A {u \D + E) n ^\a l D^ ) + (l-a 1 )E M ) ••• \v) du 
Jo du 

= E/ HA iu) ■■■ (aAD-D^ + il-a^E-E^]) •■■ \v) du 
j=i 

+ / (w\\A - A {0) ](D + E)* 11 - 1 ■■■ \v)du. (B.12) 
Jo 



On the right hand side of ( B.12 ) we regard (w\ and \v) as elements (of norm 1) of 



and note that D, E, D^°\ E<®, and D + E are bounded operators on 
of norm at most 1, and that A and A^ are bounded operators from £°° to £ 1 (Z). Now 
consider one term from the sum over j in ( B.12 ), for which the matrix product will contain 



a factor of either D — or E — E^°\ Each of these may be written as the sum of two 
terms, corresponding to the two approximations made in ( |B.2| ) or ( |B.3| ). The first of these 
terms contains a factor L — I, which we carry to the right until it reaches, and annhilates, 
\v); L — I commutes with D + E and its commutator with D^> or E^ has £°° operator 
norm which is 0(e). The second error term also has £°° operator norm which is 0(e). The 
contribution to ( |B.12|) containing A — A^ is treated similarly, and we derive an overall 



bound for ( |B.12|) which is a constant multiple of m 2 e. 
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FIGURE CAPTIONS 



Figure 1. A phase diagram for the asymptotics of the shock in the ASEP. For n >> 1, 
(r n ) ~ p + + Cn 1 exp — nn, with k = — log(l — a — 6 + ax + fe/rr) and 7 = 1 in region I and 
k — — log(l — (t/cl — V^) 2 ) and 7 = —3/2 in region II, and with C < in regions I and II a 
and C > in region lib. 

Figure 2. Profile (r n ) for p_ = 0.3, p+ = 0.6: x = 0.0 (•); x = x* = 0.2857 (x); x = 0.6 
(+). The dashed line is at height (p_ + p+)/2. 
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